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Correlation properties of the maps of the NVSS survey 
and WMAP ILC 
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In this paper we study one-dimensional sections of the maps of WMAP ILC and of the NVSS 
survey on scale lengths of 0.75, 3, 4.5, and 6.75 degrees and analyze the correlation properties of 
the sections. On these maps we identify the domains where the absolute value of the correlation 
coefficient exceeds 0.5. The catalog of such domains is presented. It is shown that the number of 
the domains agrees with the number of such domains on simulated maps and this fact may be 
indicative of just statistical agreement of the arrangement of the domains considered. 



PACS: 98.70.Dk, 98.70.Vc, 98.80.-k, 95.35.+d 



1. INTRODUCTION 

With the advent of high-precision millimiter-, 
centimeter-, and decimeter- wave radio-astronomical 
surveys it became possible to apply correlation stud- 
ies to test the self-consistentcy of the cosmological 
models. The opening of public access to one-, three-, 
and five-year accumulated data on the measurements 
of cosmic microwave background (CMB) on the entire 
sky sphere [H El Gfl H 0[7] conducted using WMAF0 
(Wilkinson Microwave Anisotropy Probe) became a 
real breakthrough in this field. To reconstruct the 
CMB signal from multifrequency observations, it is 
used the method of inner linear combinations (ILC — 
Internal Linear Combination) of background compo- 
nents [1], which yielded a map of microwave back- 
ground, also known as the ILC map, which is used 
to analyze low harmonics with multipolc numbers 
I < 100. The ILC was constructed from observations 
made in five channels: 23 GHz (the K band), 33 GHz 
(the Ka band), 41 GHz (the Q band), 61 GHz (the V 
band), and 94 GHz (the W band). 

Quite a few papers were devoted to analyzing the 
statistical properties of the signal in the reported ILC 
map and discussing the non-Gaussian nature of the 
data of this map revealed using various methods in- 
cluding phase analysis [El M 033 E] , Maxwellian mul- 
tipolc vectors PUCES], wavelets [l4l [TBI [T6l [17] . and 
Minkowski functionals [18] [19] . 

Non-Gaussian features extensively discussed by 
the researchers include the so-called Cold Spot (here- 



1 http://lambda.gsfc.nasa.gov 



after referred to as CS) with the size of about 10° 
(Galactic coordinates I = 209°, b = —57°) found in 
the distribution of cosmic microwave background in 
the Southern Hemisphere. The CS, which is statis- 
tically identified in the distribution of CMB, is one 
of the features of the map that are inconsistent with 
the hypothesis of uniform Gaussian background fluc- 
tuations. It was initially pointed out as a feature 
that deviates from Gaussian statistics when apply- 
ing wavelet analysis [T31 [201 [111 HZ] to the data of the 
first and third years of observations of the WMAP 
mission. 

A decrease of the space density of radio sources 
was found later in the NVSS maps |21| . and this 
discovery led the researches to conclude that there 
must exist a giant 140-Mpc-large void at a redshift 
of z < 1, which produces the gravitational anomaly, 
which causes the integrated Sachs-Wolfe effect [22] 
and shows up as the CS. The researchers also dis- 
cussed the possibility of explaining the spot in terms 
of the Sunyaev-Zeldovich effect [23] due to the scat- 
tering in the Eridanus supergroup of galaxies [161124] . 

Exotic interpretations of the effect included the 
texture — a topological defect during the phase tran- 
sition in the early Universe [551 US] an d Bianchi VII/j 
anisotropic cosmological model [27] . 

At the same time, serious arguments were ad- 
vanced suggesting the map contains residual contri- 
bution of Galactic background components, which 
produces the observed deviation from the Gaussian- 
ity. This contribution can show up in the form of the 
earlier discovered relation between the cleaned mi- 
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crowave background map and Galactic components 
of radiation [28l [29] and may also specify the proper- 
ties of low-order multipoles I < 20, resulting in their 
unstable reconstruction [3D]. In particular, these mul- 
tipoles may explain the specific features of the CS, 
which demonstrate the deviation of the statistics of 
clusters of signal oscillation peaks in the vicinity of 
the spot [31] — the increase of the number of positive 
peaks. 

An independent analysis of the spot's properties 
found in the region with close coordinates on maps 
and in radio-source counts in the NVSS survey [32] 
showed that the Cold Spot whose giant size and ex- 
istence are difficult to explain in terms of the cosmo- 
logical ACDM-model may be a common statistical 
deviation due to systematic effects [53] . 

In this paper we decided to verify the correla- 
tion properties of the NVSS and ILC maps on an- 
gular scales of 0.75, 3, 4.75, 6.75, and 9.75 degrees 
and to try to find and analyze CS type spots. We 
chose the angular scales to be much greater than the 
size of the smoothing window of the ILC map (1°), 
but smaller than 10° in order to preserve a sufficient 
number of pixels for statistical analysis. Our statisti- 
cal analysis also included a sub-degree angular scale 
(45' = 0.75°). In this paper we use the method of 
cross-correlation within selected sections proposed by 
Khabibullina et al. [331 131] and employed to analyze 
low-frequency spatial fluctuations of microwave back- 
ground in the domain covered by the RZF survey per- 
formed with the RATAN-600 radio telescope [3B] 137] . 
To assess the results obtained, we simulated the maps 
by generating random Gaussian fluctuations with the 
power spectrum as implied by the ACDM cosmologi- 
cal model. 

2. THE DATA 

The ILC map bears the data on microwave back- 
ground with a resolution of up to £ < 100 (about 1°). 
This map was reconstructed based on multifrequency 
observations in five different wavelength intervals 
made over the entire celestial sphere and available 
from the web-site of the WMAP mission. 

The NRAO VLA Sky Survey (NVSS) [32; is cur- 
rently the most complete high-sensitivity and large- 
area survey in the Northern sky. It was performed 
with a VLA at a frequency of 1.4 GHz from 1993 un- 
til 1996 and covers the entire Northern sky above the 
declination of 5 = -40° (33884D°or 82% of the sly 
sphere). The survey is extensively used for various 
statistical cosmological studies. The catalog contains 
a total of 1.8 x 10 6 sources and, according to its de- 
scription, is 99% and 50% complete down to the flux 
densities of 5i.4 G h z = 3.5 mJy and 2.5 mJy, respec- 



tively. The survey was made in the D-configuration 
of the VLA and the half-power size of the synthesized 
beam, which determines the resolution, was equal to 
about 45 arcsec. The data obtained as a result of the 
NVSS survey are available at the NRAO web-sit^, 
via the SkyView virtual telescope d, and from the 
CATS database0 [381 [39]. 

2.1. Properties of the NVSS catalog and the 
maps 

We analyzed the maps of the NVSS survey by con- 
structing one-dimensional sections with the smooth- 
ing due to the beam as described by Majorova [40] . 
Let us now consider two types of NVSS data: (1) the 
initial maps of the survey and (2) the catalog of ra- 
dio sources. In second case, we also set the threshold 
for elimination of sources from the scans as a free 
parameter, thereby using the possibility of forming 
scans that are free of the influence of powerful radio 
sources. 

To compare the sections constructed by differ- 
ent methods, we compute the correlation coefficients 
for: (1) scans based on the catalog and maps and 
(2) scans computed for two-dimensional maps and 
one-dimensional sections running through the center 
of the selected area. 

In first case the correlation coefficient between all 
the sections exceeds 0.99, i.e., the sections computed 
based on catalog data and those based on the maps 
can be viewed as identical for our purposes, making 
our work much easier and allowing the source lists 
from the catalog to be used for the analysis. In lat- 
ter case the correlation coefficient was 0.85 or higher, 
which we consider to be satisfactory and acceptable 
for fast computations and for the corresponding anal- 
ysis. 

To analyze the statistical correlation properties, 
we transformed the NVSS maps into a set of com- 
bined one-dimensional sections (M-scans) at given 
declinations with right ascensions ranging from h to 
24 h . We chose the declination step between the sec- 
tions to be equal to the size of the chosen pixel on 
the sphere. The values in each pixel in modified one- 
dimensional NVSS sections are equal to the mean- 
squared source fluxes within the area of given size 
centered on the given pixel: 

M * = ^rE^' P«(r)eRi, (l) 



2 http://www.cv.nrao.edu/nvss/ 

3 http://skyview.gsfc.nasa.gov 

4 http://cats.sao.ru 
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where N r is the number of sources in the pixel (ele- 
mentary sky area) of the given size of the NVSS sur- 
vey; Ri is the pixel number i of given size; S Ki is the 
flux of the radio source number Kj in pixel and 
P Ki (r) is the position of source number Ki . The pixel 
size Ri is specified by the correlation scale 0.75, 3, 4.5, 
6.75, and 9.75°. Such a characterization of the NVSS 
map makes it possible to describe the distribution of 
sources in the given sky area, which varies from pixel 
to pixel and, in particular, allows statistical verifica- 
tion of the presence of hidden fluctuations due to the 
incompleteness of the WMAP catalog [2ll4T|[42]. We 
exclude from our analysis the area in the NVSS sur- 
vey at declinations 8 < —40° and the sections that 
cross it. 

To produce sections of the ILC map we use the 
mapcut procedure of the GLESP package (43] , which 
allows sections of the map to be made at the given 
declination with the given step in right ascension. Be- 
fore performing this operation, we generated maps on 
the celestial sphere using the coefficients of spherical 
harmonics from the ai m set with the given resolu- 
tion by applying the cl2map procedure from the same 
package with the function of map synthesis: 

AT(M)=£E 

t m=l 

Yt,-m(0,(t>)) , (2) 

where Yg_ rn are the spherical functions; £ is the num- 
ber of the multipole; m is the mode of the multipole 
I; (8, 4>) are the polar coordinates, and a^ m are the 
coefficients of spherical harmonics satisfying the fol- 
lowing condition: 

Y lt _ m (9,4,) = (-l) m Y i * m (9,4 > ), 

a t , m = (-l) m <- m ■ (3) 

3. CORRELATION ANALYSIS 

The adopted approach is based on a two-steps one- 
dimensional correlation analysis between the M scans 
of the NVSS and the sections of the ILC map per- 
formed in order to search for correlation features on 
various scales. 

At the first stage we selected correlated scans of 
these two maps and identified the cases where the 
level of correlations was higher than the average level 
estimated from simulated records (Fig. 1). Although 
such a low correlation level makes the concept of cor- 
relation between the scans sounds problematic, we 
nevertheless took into account the fact that small- 
scale correlations, if present, must show up above the 
Gaussian noise and rise the correlation level K above 
the statistically expected value. 



We use standard method to compute the correla- 
tion coefficients K t for one-dimensional sections: 

_ COv(x ILC ,t, XNVSS,t) 

*M — 

&ILC,tO'NVSS,t 

n 

J2 ( x i.ILC.t — %ILC,i)( x i,NVSS,t ~ %NVSS,t) 

= ~ ; — ; ' (4) 

where XijLC,t is rth element of the one-dimensional 
section of the ILC map written in the form of vec- 
tor xiLc.t for the given scale or coordinate inter- 
val t; Xi^vss,t is the similar quantity, where vector 
x NVSS,t is based not on a section of the ILC map, but 
on the M scan of NVSS computed by formula (1); 
xihc.t and Xf g d,t are the data vectors for the sections 
of the ILC maps and M sections of NVSS respec- 
tively and <riLC,t and <JNVSS,t are the corresponding 
variances. 

The second stage consisted in analyzing the prop- 
erties of the selected sections in order to search 
for intervals with correlation coefficients |K| > 0.5. 
Let us point out from the start that if the proce- 
dure involves few points, correlations between the 
M-sections and NVSS data can always be found 
even in two random records. Our simulation of 
pure Gaussian noise using procedure fpatt of the 
FDPS reduction system [331 [35] showed that the op- 
timum number of counts such that the correlation 
coeffciient for two random records is |K| < 0.5 should 
be no less than 26. We adopted this number of counts 
as our minimum interval in pixels to search for signif- 
icant correlations (|K| > 0.5) in the records selected 
at the first stage. Figure 2 shows an example of such 
a record, which contains an interval with |K| > 0.5. 

4. RESULTS 

Due to the insufficient number of pixels in scans we 
did not use the data on for angular scales 9.75° and 
thus excluded regions with the sizes on the order of 
that of the Cold Spot. We analyze the sections and 
maps, constructed for the angular scales of 0.75, 3, 
4.5, and 6.75°. Tables 1-4 list the central coordinates 
of the intervals and the correlation coefficients. 

The number of events with K < —0.5 (34 of 53, 
19 of 34, 7 of 10, and 7 of 13 for the angular scales 
of 0.75, 3, 4.5, and 6.75° respectively) exceeds a bit 
the number of events with K > 0.5. This result could 
be explained by invoking the Sachs- Wolfe effect how- 
ever, our subsequent statistical simulation shows that 
such deviations are within the limits of Gaussian scat- 
ter. 

To test the statistical significance of the results 
with the Monte Carlo technique we used GLESP 
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Table 1: Coordinates of the centers of the areas with the correlation coefficients |K| > 0.5. The coordinates: 
right ascension a in hours and declination S in degrees (for the epoch of 2000.0) are accurate to about half the 
pixel size. The size of the pixel side is 0.75° 
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407 


-32 


87 


-0 


58 


8.387 
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Table 2: Coordinates of the centers of the areas with the correlation coefficients |K| > 0.5. The coordinates: 
right ascension a in hours and declination S in degrees (for the epoch of 2000.0) are accurate to about half the 
pixel size. The size of the pixel side is 3° 
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Table 3: Coordinates of the centers of the areas with the correlation coefficients |K| > 0.5. The coordinates: 
right ascension a in hours and declination S in degrees (for the epoch of 2000.0) are accurate to about half the 
pixel size. The size of the pixel side is 4.5° 





5° 


k 


a h 


8° 


k 


a h 


5° 


k 


21.608 


-32.12 


0.52 


15.218 


18.88 


-0.52 


22.444 


35.38 


-0.51 


12.300 


0.12 


-0.50 


16.925 


35.38 


-0.50 


23.213 


57.12 


0.57 


22.500 


0.12 


0.55 


18.764 


35.38 


-0.55 








13.950 


18.88 


-0.57 


21.708 


35.38 


-0.51 









Correlation properties of the maps of the NVSS survey and WMAP ILC 



393 



Table 4: Coordinates of the centers of the areas with the correlation coefficients |K| > 0.5. The coordinates: 
right ascension a in hours and declination d in degrees (for the epoch of 2000.0) are accurate to about half the 
pixel size. The size of the pixel side is 6.75° 
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Table 5: The number of events where the correlation coefficient exceeds K = 0.5 in the "ILC-NVSS" scans 
and in the average coefficient inferred from simulated data. The estimates are made in the scans with the same 
declination 



Angular scale, 
degrees 



0.75 
3.0 

4.5 
6.75 



Number of 
"ILC-NVSS" events 




Average number of 
"model-NVSS" events 



51.2±12.1 
22.1±6.0 
20.0±6.1 
15.3±5.4 



package [17] to generate 50 simulated random maps 
of microwave background for different angular scales 
with a power spectrum corresponding to ACDM cos- 
mology. We constructed our maps with pixelization 
on different angular scales similar to the pixelization 
of the maps studied. On each map we identified the 
scans at the declinations listed in the previous table 
and computed the number of events with the cor- 
relation coefficient |K| > 0.5. We list the computed 
number of events in Table 5. 

The results listed in Table 5 indicate that no sta- 
tistically significant deviations can be found between 
the NVSS map — ILC CMB map correlation and the 
NVSS map — random CMB maps correlation on an- 
gular scales of (1-7)°. However, we identified the do- 
mains listed in the catalog (Tables 1-4), where the 
behavior of M-scans resembles that of CMB. It is not 
improbable that the correlation between the positions 
of some spots from the catalog of celestial areas may 
be due to a real physical effect, although we would 
rather conclude that most of the correlated (anticor- 
related) spots (or perhaps even all spots) can be ex- 
plained as simple statistical coincidences. 

5. CONCLUSIONS 

We performed a correlation analysis of the CMB maps 
and of the modified NVSS map smoothed at different 
angular scales using the method of one-dimensional 
sections. We constructed the modified NVSS map 
using the distribution of the mean squared flux of 
sources within an area of given size centered on the 



pixel specified by the selected area. We used the data 
obtained as a result of our search for correlations to 
catalog the regions with |K| > 0.5 at various angular 
scales. We simulated the CMB maps using the Monte 
Carlo method and the angular power spectrum of the 
ACDM cosmology as our input parameter. We per- 
formed a similar correlation analysis for 50 generated 
CMB models. We found that the statistical properties 
of the correlations on the angular scales studied (0.75, 
3, 4.5, and 6.75°) does not differ from those of random 
maps. This result leads us to conclude that the cor- 
relations found may be due to statistical coincidence. 
However, it is of interest to perform two-dimensional 
correlation pixelization (mapping) of the sky using 
catalogs for other wavelength intervals. We plan to 
perform such an analysis in our next paper. We are 
going to, in particular, to look for the Sachs-Wolfe 
effect. 
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Figure 1: Correlation coefficients for different angular scales as a function of declination S. The correlation 
scales (0.75, 3, 4-5, and 6.75 degrees) are indicated in the figures. Gray shade shows the allowable level of cross 
correlations as computed based on the data for 50 simulated maps for KCDM cosmology. 
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